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and H. O. Kreiss 

Department of Mathematics 
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ABSTRACT 

It is shown that the dilatational terms that need to be modeled in compressible turbulence 
include not only the pressure-dilatation term but also another term - the compressible dis- 
sipation. The nature of these dilatational terms in homogeneous turbulence is explored by 
asymptotic analysis of the compressible Navier-Stokes equations. A non-dimensional param- 
eter which characterizes some compressible effects in moderate Mach number, homogeneous 
turbulence is identified. Direct numerical simulations (DNS) of isotropic, compressible tur- 
bulence are performed, and their results are found to be in agreement with the theoretical 
analysis. A model for the compressible dissipation is proposed; the model is based on the 
asymptotic analysis and the direct numerical simulations. This model is calibrated with 
reference to the DNS results regarding the influence of compressibility on the decay rate of 
isotropic turbulence. An application of the proposed model to the compressible mixing layer 
has shown that the model is able to predict the dramatically reduced growth rate of the 
compressible mixing layer. 

l'j'jijg research was supported by the National Aeronautics and Space Adminstration under NASA Con- 
tract No. NAS1-18605 while the authors were in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1 Introduction 


When the Mach number of a turbulent flow increases, the fluctuations in the thermodynamic 
variables - density, temperature and pressure - become progressively more important. Also, 
the velocity field can no longer be assumed to be solenoidal when the flow Mach number 
is significant. Turbulence modeling for compressible flow, therefore, has to account for the 
additional correlations involving both the fluctuating thermodynamic quantities and the fluc- 
tuating dilatation. In low-speed flows too, significant fluctuations in density and dilatation 
can occur in various situations, such as, the mixing of fluids with different densities, turbu- 
lent combustion, and turbulent boundary layers with strongly heated walls. This paper is 
concerned with only high-speed flows. The role of thermodynamic and dilatational fluctu- 
ations in low-speed flows is probably different from that in high-speed flows; for example, 
a supersonic shear layer at Mach 3 shows significant reduction in growth rate relative to 
the incompressible shear layer, however, a low-speed, variable-density shear layer having the 
same density difference as the Mach 3 shear layer exhibits a relatively mild change in growth 
rate with respect to its constant-density counterpart. 

Among the various additional correlations introduced into the problem due to compress- 
ibility, only the class of correlations involving the fluctuating dilatation is considered here. 
The need for modeling the pressure-dilatation is generally accepted; we show, however, that 
there is another dilational correlation - the compressible dissipation - which also merits 
attention. 

According to Morkovin’s hypothesis (Morkovin 1964, Bradshaw 1977), direct compress- 
ible effects on the turbulence may be ignored when the ratio of the root-mean-square (r.m.s.) 
density fluctuations to the mean density is small. Consequently (Bradshaw 1977), variable 
mean density extensions of incompressible turbulence models are expected to give good 
results in turbulent boundary layers with the free-stream Mach number M < 5, and in 
compressible jets with M < 1.5. Apart from the intensity of the density fluctuations, there 
is another, related, indicator of the intrinsic compressibility of high-speed turbulence, the 
turbulent Mach number M* = g/c, where q 2 is twice the turbulent kinetic energy, and c is the 
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local mean speed of sound. Using asymptotic theory and direct numerical simulations, we 
show that the compressible dissipation is naturally related to the turbulent Mach number. 

The presence of shock waves is an important feature that distinguishes the high-speed 
flows from the low-speed ones. It is known that the interaction of a shock wave with a 
turbulent boundary layer leads to significant increase in turbulence intensity and shear stress 
across the shock (Sekundoz 1974; Mateer, Brosh and Viegas 1976; Delery 1981). Some of the 
basic mechanisms underlying the shock wave/ turbulence interaction have been investigated 
through the numerical solutions (Zang, Hussaini and Bushnell 1976) of Euler equations. 
Such compressibility effects may preclude successful extension of incompressible turbulence 
models to include compressibility solely through the variability of the mean density. 

The paper is organized as follows. In Section 2 the dilatational terms that need to be 
modeled in the Reynolds stress transport equation are formally obtained. In Section 3 the 
dilatational terms are analyzed by an asymptotic theory; the main result of this section 
is the identification of a non-dimensional parameter F which is approximately equal to 
unity for low Mach number, compressible, homogeneous turbulence. In Section 4 results of 
three-dimensional direct numerical simulations (DNS) of moderate Mach number, isotropic 
turbulence are presented and shown to be in good agreement with the theoretical findings 
of Section 3. In Section 5 a model for the compressible dissipation, which is based on the 
asymptotic analysis and the DNS, is proposed; the model is calibrated witF reference to 
the DNS results on the decay rate of compressible, isotropic turbulence; and an application 
by Sarkar and Lakshmanan 2 of the new model to the compressible shear layer is briefly 
considered. Conclusions are presented in Section 6. 

2 Dilatational terms in the turbulence transport equa- 
tions 


In this section, we identify the correlations involving the fluctuating dilatation that need to 
be modeled in the Reynolds stress transport equations. It is shown that in addition to the 
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well known pressure-dilatation, an additional term, the compressible dissipation, needs to be 
modeled. 

The compressible Navier-Stokes equations, along with an equation of state, govern the 
behavior of the density p, the velocity u i} the temperature T and the pressure p in a high- 
speed, compressible flow. When the compressible flow is turbulent, an averaged form of 
the compressible Navier-Stokes equations is usually considered, wherein the instantaneous 
variables are decomposed into a mean and a fluctuating part, and the governing equations are 
averaged in order to yield equations for the mean variables. Usually Favre averages (density- 
weighted averages) are used for the velocity and temperature, while conventional Reynolds 
averages are used for the pressure and density; primarily, because such a combination leads to 
a simpler representation of the temporal derivatives and the convective terms in the averaged 
equations. We employ the above-mentioned approach too, and decompose the field variables 

as follows, 

Ui = Ui + 1i' p — P + p" 

T = T + T' p = p^p" 

The overbar denotes the conventional Reynolds average and the superscript " denotes fluc- 
tuations with respect to the Reynolds average, while the overtilde denotes the Favre average 
and the superscript ' denotes fluctuations with respect to the Favre average. The Favre 
average ^ of a field variable <f> is a density-weighted Reynolds average; 

<f> = P<f>/P 

We consider a second-order turbulence closure where in addition to the mean equations, 
transport equations are included for the Favre-averaged Reynolds stress u,u J and the turbu- 
lence dissipation rate e. The exact transport equation for u'u' is, 

2 

d^pu'u'j) + (pu k u' i u , j ) t k = Pij - Tijk,k + Ilij - ptij + ^P"u'k,k S i3 

-vJiP, j - UjP.t + u'i^jk.k + Ufrk,k (!) 
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where 


Pa 


-p(uy k u jtk + u'jU' k Ui tk ) 


= 


n a 


P e ij = 


pu[u' jU ' k + tfRs ik + F^Sik) - RRf k + ^FTk) 

^Tt + ^Kk 


In Eq. (1), Pij denotes the production, denotes the diffusive transport, II ^ denotes the 
deviatoric part of the pressure-strain correlation, and €{j denotes the turbulent dissipation 
rate tensor. The term is commonly believed to be isotropic at high turbulence Reynolds 
number; thus, for high Reynolds number turbulence, is modeled as, 

2 r 

~ 3 j 

where the turbulent dissipation rate e is defined by the expression, 


pe = a ft u ,r • 


We note that the conventional Reynolds average of the Favre fluctuation u[ is non-zero; in 
fact, u[ is related to the turbulent mass flux p /; u" by the expression 




At first glance, it appears that the only term in Eq. (l) which contains the fluctuating 
dilatation d" — u kk is the pressure-dilatation p"d". However, we show below that in a 
compressible flow there is another term containing d", which has its origin in the turbulent 
dissipation rate e. The viscous stress (J lJ in a compressible flow is given by 

°*j = p{ u i,j “b u j,i) — j 

where we have assumed that the bulk viscosity is zero. Assuming constant viscosity, the 
following expression for the turbulent dissipation rate is obtained: 


P e = 


= M(2 «--d" 2 ) 
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where the fluctuating strain rate s kl = (u kl + u" k )/2. Even if the viscosity is not assumed 
constant, standard order of magnitude estimates lead to the following expression for the tur- 
bulent dissipation rate, which is asymptotically exact for high turbulence Reynolds number: 

TK = - P*) (2) 

Let us denote the fluctuating vorticity tensor by w kl = (itjj'j — 2 and the fluctuating 

vorticity by u>". On substituting the relationship 

into Eq. (2), we obtain 

pe = p[2wl l w'l l + 2ul tl u'{ k - ^d" 2 ) (3) 

The scalar iu” k satisfies the following equation (which may be verified by inspection): 


n" 'll" 
u k t i u i,k 


= { u k u l)M ” 2 W U u i),i + 


( 4 ) 


For homogeneous turbulence, Eq. (4) becomes the following rather simple expression for 


^i" 

u k,i u i,k 


u k t k u i t i 


= d " 2 


( 5 ) 


For inhomogeneous turbulence, using standard order of magnitude estimates, Eq. (5) may be 
shown to be asymptotically correct for high turbulence Reynolds number. On substituting 
Eq. (5) into Eq. (3), we obtain 

pe = *25555+ fS*) 

= p(u/M + p" 2 ) ( 6 ) 


Thus, we have shown that for compressible turbulence the dissipation rate may be de- 
composed into two parts - the solenoidal dissipation e,, and the compressible dissipation e c , 
that is, 

pe = pe, + pe c ( 7 ) 
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where 


pe, = pwW ( 8 ) 

and 

pe c = ^JLd," 2 (9) 

Equation (7) is asymptotically exact for turbulence with high Reynolds number (which is 
of practical interest) and is exact for constant viscosity, homogeneous turbulence (which 
corresponds to the direct simulations discussed later). 

The quantities e c and e t are respectively called the compressible dissipation and the 
solenoidal dissipation because the asymptotic analysis of the next section shows that for low 
Mach number compressible turbulence e c varies on a fast compressibility time scale relative to 
e,. Only e c is directly affected by changes of compressibility indicators such as the turbulent 
Mach number while the fluctuating vorticity field and thereby e, is relatively unaffected 
by such changes in moderate Mach number turbulence. It should be noted that the direct 
numerical simulations also show that moderate compressibility affects e c and not e M . 

Zeman (1989) has also independently used a similar decomposition of the dissipation rate 
into a solenoidal and a compressible part. Zeman considers the presence of eddy shocklets 
which are assumed to augment only the compressible dissipation, bypassing the solenoidal 
energy cascade. We, on the other hand, identify the compressible and solenoidal parts of 
the turbulent dissipation by asymptotic analysis of compressible turbulence and validate the 
decomposition with direct numerical simulations. 

For polyatomic gases, the bulk viscosity may be comparable in magnitude to the shear 
viscosity p and lead to an extra dissipation which has a functional form similar to that of e c . 
The additional turbulent dissipation due to the non- negligible bulk viscosity can be easily 
modeled in the same way that the compressible dissipation e c is modeled in Section 5. 
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3 Low Mach number asymptotics 


The dimensional variables which are denoted by superscript * are non-dimensionalized as 
follows: 


p=P i 


u 


u = 


— . t = 


t*u* 


u’ 


p = 


Pr 


T = 


i; 

TL 

fl 


where 1* and u* denote a characteristic turbulence integral length scale and a turbulence 
velocity scale, respectively; /?*, p* , and T* — p*/Rp* denote reference values for the den- 
sity, thermodynamic pressure, and static temperature, respectively; and R denotes the gas 
constant. Reference values for the kinematic viscosity and the thermal diffusivity are respec- 
tively denoted by v* and a*. After using the above non-dimensionalization, the compressible 
Navier-Stokes equations take the form, 


d t p + Uip t i = -pui'i 
pd t u l -(- pujMi t j — 


1 1 


pd t T + pujTj — —(7 l)? 11 *,* 


Ml 


-qi,i + 7(7 - 1 )-Rr <7 a u i.j 


( 10 ) 

(u) 

( 12 ) 


Pr T Re T 11,1 ‘ ,w ' Re r 
The variables cr lJ and q; denote the viscous stress tensor and the heat flux, respectively. 
The non-dimensional parameters appearing in Eqs. (10)-(12) are the Mach number M r = 
u */i 'J'yP* I P *) > Reynolds number Re r = u* l* / is* and the Prandtl number Pr T = v r /cc T . 

We now consider homogeneous, compressible turbulence and adopt the approach of Er- 
lebacher et al. 2 , in which the velocity is split into an incompressible, solenoidal velocity u[ 
and a compressible velocity uf ; and the pressure change with respect to the reference pres- 
sure is correspondingly split into an incompressible pressure p 1 and a compressible pressure 
p c , as follows: 


Ui = u{ + uf 
p = 1 + p 1 + p c 


( 13 ) 
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( 14 ) 


The variables p 1 and u\ satisfy the incompressible problem, i.e., 


d t u{ + ujufj = -pj + 


— u r. 

Re r 


= 0 (15) 

The compressible problem, which is the set of equations for uf and p c , has been derived and 
discussed by Erlebacher et al. (1989) for the isentropic case. The isentropic, compressible 
problem can be investigated by an asymptotic analysis, wherein the Mach number M r is 
considered to be a small parameter, and the pressure is expanded in a power series with 
respect to M r . The leading order term in the asymptotic series for the compressible pressure 
p c is written as, 

P C = SP (16) 

where 8 — 0(M") and P — 0(1). 

We now consider the compressible problem on a sufficiently small time scale t C} which 
allows us to make the acoustic truncation of the governing equations and thus neglect the 
convective and viscous terms. (The definition of tc will be made precise later.) The problem 
for the compressible fluctuations simplifies to the following set of equations: 


s,P“ + = 0 


(17) 


a ^ y +^ p:!=0 


(18) 


The subscript 0 is used to denote the initial value 0) of a variable <f>(xi y i); for example, 


(uf)"(® i} 0) = («f") 0 (*.) 

P\*i, 0) = (P") 0 (xi) (19) 

The initial compressible velocity field (uf) 0 satisfies the conditions, 


V x .Ml") o = 0 

V(uf)o = (d") 0 
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Thus the initial value for the compressible velocity field is chosen to be irrotational and 
dilatational, while the initial value for the incompressible velocity field is chosen to be rota- 
tional and solenoidal; any arbitrary choice of the initial velocity field is amenable to such a 
Helmholtz decomposition which is unique for homogeneous flows. 

Let us denote the vorticity V X (uf )" associated with the compressible velocity by (us?) , 
and the dilatation V ' ( u f )” associated with the compressible velocity by d". On taking 
the curl of Eq. (18), and making use of the initial condition (wf )J = 0, it follows that 
(wf) W (i{, t) = 0. Thus the compressible velocity remains irrotational under the acoustic 

truncation of the governing equations. 

After some manipulation, Eqs. (17) and (18) yield, 


d tt P" 


M} 


n = o 


( 20 ) 


dud"- W d% = 0 


( 21 ) 


The above wave equations for the pressure and the dilatation are coupled through the initial 
conditions; the initial conditions for Eq. (20) are 


while those for Eq. (21) are, 


P”(*i,0) = (P")o 
d,P"(x i,0) = -j(d")o 


d"(*i,0) = (d'% 
9,d"(xi, 0) = 


( 22 ) 


(23) 


The explicit appearance of M, in Eqs. (20) and (21) is removed by resealing time through 
the transformation, 

(24) 


r = 


m t 
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After using Eq. (24) to rescale time, the equations for P" and d" take the form, 

d TT P"-P% = 0 
d TT d"-d% = 0 

The initial conditions are 

n*.,o) = (/>") o 

SrP"(Xi, 0 ) = 


(25) 

(26) 


(27) 


and 


d'\xi, 0) = (d'% 

drd'X*, 0) = -M;(P") 0 

where the quantity M* is defined by the expression, 

Mr= Wr 

We define the turbulent Mach number M t as, 

M t = q/c 

where 

q 2 = 


(28) 


(29) 


(30) 


and c is the mean speed of sound. We can rewrite Eq. (29) in terms of the turbulent Mach 
number M t as follows, 

Sq Cf 

(31) 


m;= _ 

7 M t c 


where c T = c*/u* = 1JM T is the non-dimensional, reference speed of sound. 

The system of Eqs. (25) and (26) for the pressure P" and the divergence d" are now 
solved using Fourier transforms. The Fourier transforms P(ki) and d(k t ) are defined below, 

Hk, t ) = J e tkiXi P"{xi , r) dxi 

r) = J e ik ' x ‘ d"(x t) t) dx t 
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where ki denotes the wavenumber vector. In Fourier transform space, Eqs. 
take the form, 


d rT P + k 2 P = 0 
d„d + k 2 d = 0 


(25) and (26) 

(32) 

(33) 


while the initial conditions become 


and 


P{K o) = 
d T P(ki, 0) = 


Po 

d o 


(34) 


d{ki, 0 ) = do 

d T d(k i} 0) = M*k 2 P 0 (35) 


It is a simple matter to obtain the solution of Eq. (32) and Eq. (33) which satisfies the initial 
conditions Eq. (34) and Eq. (35). The solution is 

P(ki, t ) = P 0 cos kr — s i n (^) 

d(ki, t) = d 0 cos kr + P 0 kM* sin kr (37) 

Eqs. (36) and (37) represent solutions (in Fourier space) for the compressible pressure 
fluctuation P" and the fluctuating dilatation d"; these solutions were obtained by analysis 
of the acoustic truncation of the compressible problem. Recalling that the transformed time 
coordinate r is related to the time t through Eq. (24), it is clear that the evolution of the 
compressible pressure and the dilatation from their initial values occurs in a non-dimensional 
time £ = 0(M r ). Thus for small M r , the compressible part of the problem is associated with 

a fast time scale. 

It is now possible to obtain solutions for the pressure variance P" 2 (r)> the pressure 
dilatation P"d"(r), and the dilatational variance d" 2 (r ). By definition, the pressure variance 
is related to the Fourier transform of the pressure as follows: 

W(r) = f" |P(*i, T)\ 2 iki (38) 

V J — oo 
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where V is the (sufficiently large) averaging volume, and |P(fc„r)| denotes the modulus of 
the complex variable P(^,t). Substituting Eq. (36) for P(Jfc„ r) into Eq. (38) leads to the 
following expression for the pressure variance, 


P" 2 (r) = I Ep t0 (k ) cos 2 kr dk + 


E d , 0 (k) 


Epdfi{k ) 

: sin 2 kr dk 

k 


where E^^k) denotes the initial value of the three-dimensional spectrum E^(k) of the vari- 
able <f>. 

We now introduce the concept of acoustic equilibrium value. Let <f> be a variable which 
evolves on the acoustic time scale tc- We denote the acoustic equilibrium value of a variable 
by subscript A - Then the acoustic equilibrium value (f) A of the variable <f> is the asymptotically 
stationary value that <f> attains after many acoustic time intervals. Of course, the acoustic 
equilibrium is a meaningful quantity only if the acoustic truncation of the equations apply for 
a sufficiently long time, in other words, only if the acoustic time scale t c is sufficiently smaller 
than the other time scales in the problem. The other relevant time scales in the problem of 
homogeneous turbulence are ; the turbulent time scale tp = k j e, where k is the turbulent 
kinetic energy and e is the turbulent dissipation rate; and the time scale associated with the 
mean velocity gradient = (^ij^ij) 1 ^ 2 - Since tc/tr = 0(M t ), and in usual shear-driven 
turbulent flows = 0(1), the acoustic equilibrium is formally realizable when M t << 1. 

It should be noted that the acoustic equilibrium value corresponds to a quasi-equilibrium 
phase during the evolution of the variable; the variable remains stationary over the time 
interval tc « t « tp. 

Mathematically, the acoustic equilibrium value of P" 2 is obtained by evaluating Eq. (39) 
in the the limit r -* oo. The Riemann-Lebesgue theorem, which states that f(k)e ikt dk — ► 0 
as < — ► oo provided / a 6 \f(k)\ dk exists, is used to help evaluate this limit, and the following 
expression for the acoustic equilibrium (P" 2 ) A of the pressure variance is obtained: 
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Let us denote the compressible portion of q 2 (which is twice the turbulent kinetic energy) 
by by definition, we have the relation, 


<& = MW)" 

The three-dimensional spectrum Ejjt) of the dilatation and the three-dimensional spectrum 
E q a (&) of the compressible portion of q 2 are related by 

Ej(k) = k 2 E${k) (41) 


Using this relation, Eq. (40) yields 


(P" 2 )* = i + ft 1 


(42) 


where F 0 denotes the initial value of the non-dimensional parameter F which is defined as, 


F = — — (43) 

MfP" 2 

An expression for the acoustic equilibrium of the compressible turbulent kinetic energy 
(?c)a is now obtained. After rescaling time by 

t 


T = 


M r 


Eqs. (17) and (18) become 


d,p"+¥g-(W = 0 


(44) 

Multiplying Eq.(44) by {2S 2 P")/(^M 2 ), multiplying Eq.(45) by 2(uf)", and adding the two 
resulting equations gives 


M* 2 d r P" 2 + dr(uf)" 2 + 2 M;[P"(u?)"],i = 0 


nf..c\in 


(46) 


where M; = «/(qM r ). Averaging Eq. (46) gives the following result for homogeneous tur- 
bulence, 

d T [q 2 c + M' T 2 W) = 0 (47) 
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Thus, the quantity in square brackets in Eq. (47), which physically represents the full non- 
dimensional turbulent energy, does not change on the acoustic time scale, and consequently 
we have 

{q 2 c) A + M; 2 (P^) a = (<£) o + M; 2 (p^) 0 (48) 

After substituting Eq. (42) into Eq. (48), we obtain 

(<1c)a = gfacOof 1 ] (49) 

where F is defined by Eq. (43). 

On dividing Eq. (49) by Eq. (42), we obtain the interesting result that the acoustic 
equilibrium of the non-dimensional parameter F is unity; 


f a - 1 (50) 

The physical significance of F is better understood by reverting to dimensional quantities 
(denoted by superscript *). After some manipulation, F may be written as 

F=- ^ f5n 

0 Pc*") 2 hP'r m 

The numerator of Eq. (51) is twice the kinetic energy of the compressible turbulence, and the 
denominator is twice the potential energy of the compressible turbulence. Thus, the result 
F a — 1 implies that at acoustic equilibrium there is an equipartition between the kinetic 
and potential components of the compressible energy. Since Eq. (50) is a consequence of 
processes occurring on the acoustic time scale tc, as long as the other time scales of the 
problem (such as kj e) are larger than t c , we have F ~ 1 for later time. Therefore low Mach 
number, homogeneous, compressible turbulence has an equilibrium structure characterized 

by 

F - 1 (52) 

We will now derive an alternative expression for F in terms of the turbulent Mach number 
M t which will be useful later. On substituting the expression for M* from Eq. (31) into 
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Eq. (43); recognizing that c and p are approximately constant on the acoustic time scale, 
and are respectively equal to their initial, references value <v and ;v; we obtain 

l 2 Mix 


F = 




(53) 


Here x denotes the ratio of compressible kinetic energy to the total turbulent kinetic energy, 
that is, 


9c 

x = 7 


and p c is the non-dimensional ratio of the r.m.s pressure fluctuations to the mean pressure, 

(54) 

P,= — 


An expression is now sought for the equilibrium value of the pressure-dilatation P"i". 
Starting with Eqs. (36) and (37), the pressure-dilatation can be related to the initial condi- 

tions of the turbulence through, 

Wd" = E Pd ,o(k) cos 2kr dk-^jj Q ^^-sm2kr dk 

4 . Ml. f kE Pi0 (k)sm2kTdk 

2 Jo 

Again using the Riemann-Lebesgue theorem to evaluate the right hand side of Eq. (55) in 

the limit r -> oo, gives (P"d") A = 0. ..... 

The acoustic equilibrium value of the dilatational variance d is obtained in a si 

manner. The expression for d" is 

d^(r) = r EaM cos 2 kr dk + M r * 2 jf" k 2 E Pfi (k) sin 2 kr dk 

+ Ml f°° kE P d,o{k) sin 2kr dk 
Jo 

while the expression for (d ,,2 )x is 




(56) 
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Thus, of the two dilatational correlations P"d" and d"\ the acoustic equilibrium of W is 
zero, while the acoustic equilibrium of the positive definite quantity d is non-zero. 

To summarize, in this section we have identified certain variables of compressible, ho- 
mogeneous turbulence which evolve from arbitrary initial conditions on a fast time scale t C) 
where i c = 0(M t k/e). The time evolution of these variables has a quasi-equilibrium phase 
m which the variable has a stationary value which we call the acoustic equilibrium value. 
We have also shown that these variables can be combined into a non-dimensional parameter 

F which, after starting from an arbitrary initial value, maintains a value of approximately 
unity; thus 


fc ' ' 

Where M ‘ = “ lhe ‘uAulent Mach number. * i. the ratio of compressible turbulent kinetic 

energy to the full turbulent kinetic energy, 7 is the ratio of specific heats, and p c = jp 

is the ratio of the r.m.s. compressible pressure to the mean pressure. 

4 Direct numerical simulation of compressible, isotropic 
turbulence 

Three-dimensional direct numerical simulations (DNS) of compressible, isotropic turbulence 
were performed on a 64 grid for a variety of initial conditions. Details of the algorithm and 
numencs are provided in Erlebacher et al. (1987). The simulations correspond to a nominal 
turbulence Reynolds number Re x (based on the Taylor microscale) of 15. 

The behavior of the compressible dissipation and the pressure-dilatation for a case with 
initial turbulent Mach number M t>0 = 0.5 is illustrated in Figs. 1-2. The compressible 
dissipation reaches its acoustic equilibrium value after a time boundary layer of 0(kM t /e), 
and then decays with a small superimposed acoustic modulation. The pressure-dilatation, 
which can be of either sign, shows a significant acoustic modulation, and tends to be more 
positive than negative. The direct simulations indicate that in the case of decaying isotropic 
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turbulence, the pressure-dilation, when averaged over the acoustic oscillations, is positive 
and smaller than the compressible dissipation. 

The asymptotic analysis of the previous section predicted that the nondimensional pa- 
rameter F = {l 2 M?x)/pl should be approximately equal to 1. The DNS show that, after 
starting from a variety of initial values, F indeed reaches a value of unity. Fig. 3a shows the 
early-time behavior of F for three representative cases; F attains its acoustic equilibrium 
value of unity in a non-dimensional time of 0{M t ). Fig. 3b, which depicts the late-time be- 
havior of F, shows that F exhibits relatively small excursions from its acoustic equilibrium 
value of unity Even though the individual quantities such as M? decrease by about a factor 
of 3 in the time interval 0.4 < (e a ) 0 t/k 0 < 2.0, the quantity F deviates from its theoretically 
predicted value of unity by less than 10%. If F is averaged over a few of its oscillations, the 
deviation of this averaged value from unity would be much smaller than 10%. 

Fig. 4 shows the DNS results on the decay of the turbulence kinetic energy k for three 
values of the initial turbulent Mach number Mt,o- The initial value of the nondimensional 
r.m.s pressure fluctuation was chosen as M* 0 and in order to eliminate the initial time 
boundary layer F was set equal to unity It is clear that an increase in the compressibility 
level tends to increase the decay rate of k. Evidently, compressibility leads to an additional 
source of dissipation for the turbulent kinetic energy. 

5 Modeling of the dilatational terms 

We will now develop models for the two dilatational terms - the pressure-dilatation, and 
the compressible dissipation - that appear in the Reynolds stress transport equations. The 
theoretical analysis indicated that low Mach number homogeneous turbulence is character- 
ized by the relation F ~ 1, and the DNS showed that F ~ 1 for turbulent Mach numbers 
at least up to M t = 0.5 (which, in free shear flows, corresponds to the mean Mach num- 
ber M being as large as 10). We will now make use of the result F ~ 1 for developing a 
simple algebraic model of the compressible dissipation e c . The model essentially relates the 
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turbulent Mach number, which is perhaps the most important quantity characterizing the 
intrinsic compressibility of high-speed turbulence, to the compressible dissipation. 

The compressible fraction of the dissipation rate % e = e c / 6 satisfies the following equation. 

y = ^ (58) 

ydiere the compressible Taylor microscale A c , the solenoidal Taylor microscale A # , and the 
compressible fraction of turbulent kinetic energy x, are defined as follows, 


A c = (qh/^ 2 (59) 

A. = (g|/^f) 1/2 ( 6 °) 

x = <&!<? ( 61 ) 


Using Eq. (58) we obtain the following expression for f €} the ratio of compressible dissi- 
pation to the solenoidal dissipation: 

t - h- 4A,2% (62) 

J '~ e.~ 3A c 2 (1 - x) 

We assume that for compressible turbulence, A c /A, = 0(1), and from Eq. (62) obtain the 
asymptotic representation for small x : 

/e = Ax + °(x 2 ) ( 63 ) 


where = 0(1). On using Eq. (57), and recognizing that p c = 0(M 4 2 ), we obtain the 
following expression from Eq. (63): 

e c = e.^M 2 + 0(M t 4 )] (64) 

where a x = 0(1). We now propose the following algebraic model for e c , which is motivated 
by Eq. (64): 

e c = c*iM t 2 e, (65) 

where cci is a constant of 0(1), whose numerical value remains to be evaluated. We note that 
a natural extension of the model Eq. (65) for large M t is to add a term proportional to M * 
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in Eq. (65). For now, we will limit ourselves to the simpler model Eq. (65). The solenoidal 
dissipation rate e, is calculated using the standard form of the incompressible dissipation 
rate transport equation. 

We also need a model for the pressure-dilatation p"d". The asymptotic analysis predicts 
that the acoustic equilibrium of the compressible pressure-dilatation P"d" is zero, while that 
of the compressible dissipation e c is non-zero. The DNS indicate that in the case of isotropic 
turbulence, except for the initial time boundary layer, the average of the pressure-dilatation 
pf'd" over its oscillations is smaller than the compressible dissipation e c . Therefore, for the 
purpose of turbulence modeling, we absorb the effect of p"d" in the model of e c . It should be 
noted that p"d" may be non-negligible if eddy shocklets are present; a model of the pressure- 
dilatation which is applicable to the eddy shocklet regime has been formulated by Lumley 
(1989). 

The model constant a : is evaluated by considering the compressible, iso-decay problem. 
After introducing the models for the compressible dissipation and the pressure-dilatation, 
and using the standard dissipation transport equation, the governing equations become 


UbfJ 

dt 

0 

(66) 

II 

— e,(l + eiiMl) 

(67) 

dk 

dt 

-e.(l + <*iMl) 

(68) 

d{e.) 

dt 

c e ’ 
~ Ct2 J 

(69) 


Since Ml = ? 2 /(c) 2 , Eqs. (67) and (68) can be combined into the following equation for A/ 2 : 

= ~jMl(l + ai Mf)[ 1 + 0.57(7 - 1 )Ml] (70) 

The term, — e c + p"d", which is the extra compressible term on the right-hand-side of the 
exact turbulent kinetic energy equation, is replaced by the model, — ajM t 2 e 4 , in Eq. (68). 

Eqs. (68), (69) and (70) were integrated with a fourth-order, Runge-Kutta scheme using 
various values for both the model constant and the initial Mach number M ti0 . The results 
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of these computations were then compared with the DNS. The model coefficient C e2 was 
chosen to be 1.83 so as to reproduce the experimentally observed decay rate in physical ex- 
periments on high Reynolds number incompressible turbulence. Since the Reynolds number 
of the simulations is somewhat low, the turbulence decays faster in the simulations relative to 
the high Reynolds number experiments. Therefore, when comparing model results with the 
DNS, rather than looking for agreement between the actual value of the decay rate obtained 
with the model and that obtained with the DNS, we look for agreement regarding the effect 
of compressibility on the turbulence decay rate. After comparing Fig. 4 and Fig. 5 it is clear 
that, as far as the influence of the turbulent Mach number on the decay rate of the turbulent 
kinetic energy is concerned, the choice of ai = 1 gives good qualitative agreement between 
the the results of the model and the DNS. Thus, the model for the compressible dissipation 

becomes 

e c — a\M*e, (71) 


where the model constant c*i = 1. 

It should be noted that the model (71) for the compressible dissipation has been applied 
by Sarkar and Lakshmanan 3 to the compressible shear layer, within the framework of a Favre- 
averaged Reynolds stress closure. Details regarding the other modeling assumptions in the 
closure and the numerical implementation of the second-order closure; and results for various 
configurations of the compressible shear layer are provided by Sarkar and Lakshmanan. Fig. 
6 is a schematic of the particular configuration of the shear layer, a few of whose results 
are given here. A high-speed stream with velocity U\ mixes with another stream with lower 
velocity U 2 . The free-stream values of the pressure, density and temperature are equal in 
the two streams. The normalized spreading rate C s defined by 


3 ICASE report in preparation 


&.V1 + U2 

6 dx ^U\ — U 2 
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is the primary variable of interest. The shear layer thickness 6(x) is defined to be the distance 
between the two points of the mean velocity profile where the mean velocity is respectively 
U 2 + 0.1 (Ui - U 2 ) and U 2 + 0.9 {U x - U 2 ). 

Fig. 7 shows model predictions and experimental data on the influence of mean compress- 
ibility on the spreading rate of the mixing layer. The mean compressibility of the compress- 
ible shear layer is characterized by the relative Mach number Mr = 2{U\ — U 2 )/(ci + c 3 ), 
where Ci and c 2 denote the free-stream speed of sound in the two incident streams. It should 
be noted that, for the shear layer of Fig. 6, the relative Mach number Mr is twice the 
convective Mach number M c as defined in Papamoschou and Roshko (1988). In Fig. 7, we 
plot the non-dimensional spreading rate Cg/(Cg) o, where Cg is the spreading rate of of the 
mixing layer and (Cg ) o is the spreading rate of the incompressible mixing layer. Though 
there is a systematic difference between the data of Papamoschou and Roshko (1988) and 
the data of the Langley curve (see Kline et al. 1981), which is a consensus representation 
of various experimental investigations, it is clear that first, the spreading rate decreases sig- 
nificantly when the relative Mach number increases; and second, after the initial decrease, 
the spreading rate is relatively insensitive to further increases in the relative Mach number. 
The prediction of the second-order closure, with the model for the compressible dissipation 
included, is in agreement with both the aforementioned trends exhibited by the experimental 
data. However, excluding the model of the compressible dissipation from the second-order 
closure leads to only a mild decrease of spreading rate with increasing Mach number. 

6 Conclusions 

Asymptotic analysis of the compressible Navier-Stokes equations has led to the identification 
of a non-dimensional parameter F which characterizes some compressible effects in high- 
speed turbulence. The variable F evolves from arbitrary initial values on a non-dimensional 
time scale of 0(M t ), attains an equilibrium value of unity, and remains approximately equal 
to unity for later time. The result F ~ 1 , which is formally valid only for turbulent Mach 
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number M t « 1, has been shown to hold in the direct numerical simulations (DNS) of 
isotropic turbulence where M t was varied between 0.01 and 0.5. 

It was established that there is another dilatational correlation - the compressible dissi- 
pation - which needs to be modeled in addition to the well-known pressure-dilatation. Both 
the theoretical analysis and the direct simulations suggest that the compressible dissipation 
is larger than the pressure-dilatation in low Mach number, homogeneous turbulence. A sim- 
ple algebraic model, which is based on asymptotic analysis and DNS, has been proposed 
for the compressible dissipation. The model, which has been applied to the calculation of 
a high-speed shear layer, was able to capture the dramatically reduced growth rate of the 
high-speed shear layer. 

The present turbulence closure, where dilatational effects are included through a simple 
model having an algebraic dependence on the turbulent Mach number, will be extended in 
the future to include transport equations for the thermodynamic turbulent statistics such as 
the density variance. The consequence of higher order extensions of the asymptotic theory 
to compressible turbulence modeling will also be explored. 
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Figure 6. Schematic of the compressible shear layer. 
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Figure 7. Application of the model for compressible dissipation to the compressible shear 
layer; from Sarkar and Lakshmanan (1989). 
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